Compositional (version 5.5)

Contour plot of the alpha-folded model in S^2: Contour plot of the $\alpha$-folded model in $S^2$

Description

Contour plot of the $\alpha$-folded model in $S^2$.

Usage

folded.contour(mu, su, p, a, n = 100, x = NULL, cont.line = FALSE)

Arguments

mu

The mean vector of the folded model.

su

The covariance matrix of the folded model.

p

The probability inside the simplex of the folded model.

a

The value of a for the $\alpha$-transformation.

n

The number of grid points to consider over which the density is calculated.

x

This is either NULL (no data) or contains a 3 column matrix with compositional data.

cont.line

Do you want the contour lines to appear? If yes, set this TRUE.

Value

The contour plot of the folded model appears.

Details

The $\alpha$-transformation is applied to the compositional data and then for a grid of points within the 2-dimensional simplex the folded model's density is calculated and the contours are plotted.

References

Tsagris M. and Stewart C. (2020). A folded model for compositional data analysis. Australian and New Zealand Journal of Statistics, 62(2): 249-277. https://arxiv.org/pdf/1802.07330.pdf

See Also

alfa.contour, compnorm.contour, diri.contour, mix.compnorm.contour, bivt.contour, skewnorm.contour

Examples

# NOT RUN {
x <- as.matrix(iris[, 1:3])
x <- x / rowSums(x)
a <- a.est(x)$best
mod <- alpha.mle(x, a)
folded.contour(mod$mu, mod$su, mod$p, a)
# }